O 
O 

(N 



We have developed software (SearchCal) that builds 
an evolutive catalog of stars suitable as calibrators within 
any given user-defined angular distance and magnitude 
around the scientific target. We present the first version 
of SearchCal dedicated to the bright-object case (V<10; 
K<5). Star catalogs available at the CDS are consulted 
via web requests. They provide all the useful information 
for selecting of calibrators. Missing photometries are com- 
puted with an accuracy of 0.1 mag and the missing angular 
diameters are calculated with a precision better than 10%. 
For each star the squared visibility is computed by tak- 
ing the wavelength and the maximum baseline of the fore- 
seen observation into account. SearchCal is integrated into 
ASPRO, the interferometric observing preparation soft- 
ware developed by the JMMC, available at the address: 
http: / /mariotti.fr 

Key words. Techniques: high angular resolution 
Techniques:interferometric- Stars: fundamentals parame- 
ters - Catalogs - Astronomical database: miscellaneous 
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Abstract. In long baseline interferometry, the raw fringe contrast must be calibrated to obtain the true visibility 
and then those observables that can be interpreted in terms of astrophysical parameters. The selection of suitable 
calibration stars is crucial for obtaining the ultimate precision of interferometric instruments like the VLTI. 
Potential calibrators must have spectro-photometric properties and a sky location close to those of the scientific 
target. 



1. Introduction 

A long baseline optical interferometer measures the spatio- 
temporal coherence (i.e. the visibility) of the target. This 
measure is directly related to the Fourier transform of the 
object's intensity map at the spatial frequencies ^, B 
being the baseline vector between two telescopes. Optical 
interferometry provides a powerful tool for determining 
the morphology of astronomical sources at high angular 
resolution. The modulus and the phase of the visibility are 
derived respectively from the contrast and the position of 
the fringes resulting from the recombination process. 

The atmospheric turbulence and instrumental insta- 
bilities induce long-term and short-term drifts that dis- 
tort the phase and decrease the amplitude of the target's 
visibility Vtarget by a factor T called the instrumental re- 
sponse. The observed visibility u, target could then be writ- 
ten as 

A* target = V target^ ■ (1) 

In order to take these effects into account and to con- 
vert the observed fringe contrast into true visibility, the 
observation of the scientific target is usually bracketed by 
observations of calibration stars. In practice, observing a 
calibration star whose visibility V ca i can be accurately de- 
duced from direct or indirect determination of its angular 
diameter leads to a determination of T: 

-p _ \Hcal | 

\vu {2) 



It is then possible to calculate the accuracy on the target's 
visibility: 

AV t l rget _ AV 2 d A M 2 1 1 
V? ~ V 2 V 2 ^V 2 V 2 ' [ ' 

'target y cal cal "target 

where Afj 2 ~ An 2 target ~ Afi 2 cal is the uncertainty on 
the measurement of the visibility amplitudes. Equation 3 
shows that the expected accuracy on the target visibility 
strongly depends on the accuracy of the calibrator visibil- 
ity. 

A calibrator is a star for which the visibility is known 
(or can be predicted) with high accuracy. It should have 
physical properties (magnitude, spectral- type, colors) and 
a sky location close to those of the scientific target, so that 
the instrumental response during the calibrator-target- 
calibrator sequence could be considered as independent 
of the object. 

The selection of suitable calibration stars is crucial 
to obtain the ultimate precision of the interferometric 
instruments. Until now, each interferometric group has 
had its own strategy for calibrating the observations ei- 
ther by using reference stars chosen case by case or 
using specific tools of selection. In 2002, Borde et al. 
published a catalog of 374 reference stars selected from 
the initial list of Cohen's spectro-photometric calibrators 
(Cohen et al. [1999) using selection criteria adapted to in- 
frared interferometry up to a 200 m baseline. More re- 
cently, an observing program (Percheron et al. 2003) was 
set up to create a list of reference stars with accurate 
measured angular diameters suitable for calibrating the 
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infrared intcrfcromctric observations of the Very Large 
Telescope Interferometer (VLTI) instruments (VINCI, 
MIDI, AMBER). To prepare interferometric observations 
with the Palomar Testbed Interferometer (PTI) and Keck 
Interferometer (KI), an interferometric observation plan- 
ning software GetCal has been developed (Boden [2003), 
including a tool to compute the visibility of potential 
reference stars taken in the Hipparcos catalog and ex- 
tracting astronomical and spectro-photometric parame- 
ters from the Simbad database at the Centre de Donnees 
Astronomiques de Strasbourg (CDS)(Genova et al., 2000). 

With the startup of long-baseline and large-aperture 
optical interferometers, such as VLTI, KI, or the CHARA 
array, and with the increase in the accuracy or in the 
range in sensitivity and angular resolution, the calibration 
of interferometric data requires definiting new strategies 
for seeking suitable calibrators and developing of selection 
tools usable by a larger community of astronomers. 

In Sect. 2, we present our method for creating a dy- 
namical list of stars fulfilling the requirements of interfer- 
ometric calibrators for a bright scientific target^ < 5). 
Section 3 briefly depicts the different sccnarii of request to 
the CDS database, in order to extract the useful parame- 
ters from stellar catalogs and to sort out the initial list of 
possible calibration stars. Section 4 deals with the major 
steps in the calculations (interstellar absorption, angular 
diameters, visibility) for each star on the list. Some techni- 
cal aspects are mentioned in Sect. 5. Finally, the current 
limitations of SearchCal and its evolution to in case of 
faint targets are discussed in the last section. 

2. Our original method 

The design of a search calibrator tool available in ASPRO 
(Duvert et al., 120021 Duchene et al., I2004[> was guided 
by the goal of creating a dynamical catalog of calibration 
stars suitable for each scientific target. The goal was to 
provide a list of potential calibration stars for which the 
visibilities are calculated from their angular diameters and 
the maximum spatial frequency (^) of the interferometric 
observation. The search for calibrators must work as well 
for long baseline interferometric observations carried out 
in the visible (V band), the near infrared (J, H, or K 
bands), or the mid infrared (N band). 

The " Virtual Observatory" techniques were adopted to 
extract the required astronomical information from a set 
of stellar catalogs available at the CDS. Compared to the 
static or closed-list approach, the merit of this strategy is 
first to take into account any enrichment of the catalogs 
by new observational data and secondly to be much more 
adapted to the limits in magnitude of the coming inter- 
ferometric facilities (VLTI with four 8 to or KI with two 
10 to telescopes). 

To minimize the effects of temporal and spatial varia- 
tions of the seeing on the calibration process, a calibrator 
must be as close as possible to the scientific target. The 
field size on the sky is defined by the maximum differ- 
ence in right ascension and declination. To be observable 



with the same instrumental configuration, the magnitude 
of the calibrator must be in a narrow range of value around 
the target magnitude in the observing photometric band. 
In order to select stars as potential calibrators, a certain 
number of astronomical parameters must be known for 
each star. These parameters are given in Tabled 

Table 1. Astronomical parameters for calibration stars 

Identifiers HIP, HD, DM numbers 

Astrometry coordinates (RAJ2000, DEJ2000), 

proper motion, parallax, 

galactic coordinates 
Spectral Type temperature and luminosity class 

Photometry magnitudes U, B, V, R, I, J, H, K, L, M, N 

Angular diameter measured or computed angular diameter 
Miscellaneous variability and multiplicity flags, 

radial velocity, rotational velocity 



An on-line interface with the VizieR data base 
(Ochsenbein et al.. I2000|) at the CDS was created to ex- 
tract astrometric and spectro-photometric parameters of 
the sources in the defined box and to obtain the initial 
list of stars (see details in the next section). This list is 
enriched by the stars present in the Catalogue of calibra- 
tors for long baseline stellar interferometry (Borde et al., 
2002) and the Catalog of bright calibrator stars for 200- 
m baseline near-infrared stellar interferometry (Merand 
et al., 2005). If available, the measured angular diame- 
ter is obtained through the data of the Catalog of High 
Angular Resolution Measurements (Richichi, Pcrcheron 
and Khristoforova. l2T)05f) . 

For each star on the initial list, calculations are made 
to correct the interstellar absorption and to compute miss- 
ing magnitudes. The photometric angular diameter and its 
associated accuracy are estimated using a surface bright- 
ness method based on the (B — V), (V — R) and (V — K) 
color index. Then, the expected visibility and its error are 
computed. 

The list of possible calibrators is finally proposed to 
the user and the final choice can be made by changing the 
selection criteria: accuracy on the calibrator visibility, size 
of the field, magnitude range, spectral type and luminosity 
class, variability and multiplicity flags. 



3. The CDS interrogation 

To build a dynamical list of stars, we chose to extract the 
information from catalogs available at the CDS. Different 
scenarii were implemented depending on the photometric 
band selected for the interferometric observations, i.e. in 
the visible (V band) or the near infrared (K band). For 
each star the astronomical parameters were extracted from 
the following catalogs: 

— 1/280: All-sky Compiled Catalog of 2.5 million stars 
(Kharchenko, 2001) 
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- II/7A: UBVRIJKLMNH Photoelectric Catalog (Morel 
et al., 1978) 

- 11/225: Catalog of Infrared Observations, Edition 5 
(Gezari et al, 1999) 

- II/246/out: The 2MASS all-sky survey Catalog of 
Point Sources (Cutri et al., 2003) 

- J/A+A/413/1037: catalog J-K DENIS photometry of 
bright southern stars ((S. Kimeswenger et al.,2004) 

- I/196/main: Hipparcos Input Catalog, Version 2 
(Turon et al., 1993) 

- V/50: Bright Star Catalog, 5th Revised Ed. (Hofflcit 
et al., 1991) 

- V/36B: Supplement to the Bright Star Catalog 
(Hoffleit et al., 1983) 

- J/A+A/393/183: Catalog of calibrator stars for LBSI 
(Borde et al., 2002) 

- J/A+A/433/1155: Calibrator stars for 200-m baseline 
interferometry (Merand et al., 2005) 

- J/A+A/431/773/charm2: Catalog of High Angular 
Resolution Measurements (Richichi et al., 2005) 

To define of the extracted data and the limitation of 
the number of returns, we used for each catalog the data 
fields defined by UCDs (Unified Content Descriptors) and 
labels and the limits on the data's values. Our strategy 
is based on two sequences of requests on the VizieR data 
base. 

3.1. Primary request 

An on-line interface with the CDS has been created to ob- 
tain the initial list of stars present in the calibrator field 
and with magnitudes according to the specified magnitude 
range. This request is done on catalog(s) called "primary 
catalog(s)" depending on the scenario. The primary cata- 
logs were selected with respect to the quality of the equa- 
torial coordinates and of the available photometry. In the 
case of the "visible" scenario, the choice was thus made 
on the compiled catalog 1/280 because of the need to have 
a reliable value for the magnitude V and precise coordi- 
nates for stars brighter than typically Vmag < 10. For the 
"near infrared" scenario, it was mandatory to have the K 
magnitude of a star brighter than typically Kmag < 5, 
and the choice was made to take the compiled catalogs 
1/225, II/7A, and 11/246 as primary catalogs. The output 
of this first sequence of requests is a list of star coordi- 
nates having magnitude values as specified in the defined 
calibrator field. 

3.2. Secondary request 

The second sequence of requests is done on the stars con- 
tained in the previous list and with the goal of extract- 
ing astrometric and spectro-photomctric parameters. The 
secondary catalogs were selected because of the relevance 
and the reliability of the parameters of interest for our 
purpose. In the current version of SearchCal, the identi- 
fiers, the equatorial coordinates, the proper motions, the 



parallaxes, spectral type, and variability or multiplicity 
flags are extracted from the 1/280 catalog. For bright 
stars, the visible photometry comes from 1/280, whereas 
infrared photometry is taken from the II/7A, 11/225 and 
11/246 catalogs. The galactic coordinates are taken from 
the 1/196 or 11/246 catalogs. The radial velocity and ro- 
tational velocity are extracted from the 1/196, V/50 or 
V26B catalogs, respectively. 

3.3. Setting the list of possible calibrators 

We then parse and merge all the results in a single ar- 
ray of stars with all the astronomical parameters. The 
catalogs are first linked according to the HD number if 
provided and with the equatorial coordinates found in the 
different catalogs if they arc coherent at the level of 1 arc 
second. The V magnitude (for V band) or the K magni- 
tude (for infrared bands) is also used to confirm that the 
star present in the different catalogs is the same. The final 
result is a single list containing stars for which the suit- 
able astronomical parameters have been extracted from 
the selected catalogs. 

4. The central engine of the calibrator's 
parameters calculation 

For the stars contained in the final list of the CDS re- 
quests, we need to compute their apparent diameters (ex- 
cept if they have already been measured) to determine 
their visibilities in the interferometric configuration. This 
is done in several stages: first we correct the photometry 
from interstellar absorption, then we compute the possible 
missing photometric data, and finally the apparent diam- 
eter is obtained from surface brightness relation. 

4.1. Interstellar absorption 

We must correct the photometric data for the wavelength- 
dependent effects of the galactic interstellar extinction. In 
the current version of SearchCal, all the calibrators are 
bright enough to have a measurement of their trigono- 
metric parallax. For each star, the visual absorption Ay 
is computed as a function of the galactic coordinates 
(longitude I and latitude b) and of the distance d. As 
all calibration stars currently selected by SearchCal have 
d < lOOOpc, we have used the analytic expression for the 
interstellar extinction in the solar neighborhood given by 
Chen et al. iTHMIl . 

The observed magnitudes mag [A] are then corrected 
for interstellar absorption using: 

mag[X]o = mag [A] - A x (4) 
Rx = A X /E(B - V) (5) 
A x = AyRx/Rv (6) 
with Ry = 3.10 and the values of Rx given by Fitzpatrick 
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4.2. Rebuilding missing photometries 

The knowledge of all the BVRIJHK photometry of the 
calibrator is useful for the interferometric observation. In 
particular, it could be important to have a calibrator with 
similar brightness to the scientific target at the observ- 
ing wavelength, and the photometry will be mandatory 
for computing the angular diameter of the calibrator. In 
addition, the value of the magnitude of the calibrators in 
some photometric bands must be known for the needs of 
some housekeeping operations (fringe tracking for exam- 
pie). 

If some photometric values are missing, we have chosen 
to compute them from "spectral type - luminosity class - 
color" relations and from existing magnitudes. Such rela- 
tions exist in the literature, but in general each of them 
covers only certain classes of luminosity, a range of spec- 
tral type, or a limited number of photometric index. To 
obtain a relation linking all spectral types of all luminosity 
classes to BVRIJHKLM photometry, we compiled the 
works of Bessel (|1979fl . Bessel and Brett l|1988j) . Fitzgerald 
fHJ70J> Johnson ifTMSjl . Leggett l(T^|l Schmidt-Kapler et 
al., lfT3g2|) . The ijTMijl and Wegner <|TOMjl . We adopted 
the Johnson photometric system according to our main 
source of accurate photometry (Morel & Magnenat ll978|l . 
The relation of Bessel (jl98rjfl , Glass H1975[l and Bessel and 
Brett ( 1988) are used to transform the other photometric 
systems in the Johnson one. Tables lA~TllA.2l and lA.3l list 
the adopted value of our "spectral type - luminosity class 
- color" relations for dwarfs, giants and supergiants stars, 
respectively. In Fig. ^ we show an example of our relation 
for the (B - V) and (V - R) color index. The [B - V) re- 
lation is also used to check the consistency of the spectral 
type extracted at the CDS. 

To check the accuracy of the rebuilt photometry, we 
compare the colors of our tables with those of stars in the 
catalog of Ducati (2002) also in Johnson filters. Figure |21 
shows an example of the O-C (difference between the color 
measured and computed with our relation) as a function 
of the spectral type. The Ducati (2002) photometry is not 
corrected for interstellar absorption, so a part of the dis- 
persion is due to reddening, which is visible for the bright 
(and then distant) OB stars. The dispersion of the O-C is 
then a superior limit of the accuracy of our computed pho- 
tometry Using only the stars in the spectral type range A 
to K seems a good compromise for estimating our accu- 
racy, since they are closer than OB stars so less reddened 
and enough bright to reduce the observational errors. 

In the case of near infrared observations, we impose 
the knowledge of the photometry in V and K, and our 
computed complementary photometry has an accuracy of 
0.1 mag or better. For visible observations, only the B and 
V magnitudes are mandatory, so that the determination 
of the missing J, H , and K has an accuracy of 0.2 mag. 




1 2 

V-R 

Fig. 1. Example of our "spectral type - luminosity class - 
color" relations used to compute missed photometry. The 
relation for dwarfs, giants and supergiants are shown in 
dark, red, and blue, respectively . Our sequences are plot- 
ted in solid lines, when the relation from which they are 
extracted are in different dashed lines. 

4.3. Determination of the angular diameter 

Since the goal was to calculate the visibility for each pos- 
sible calibrator, it was necessary to know the value of its 
angular diameter. In some particular cases this parameter 
is either a measured value taken from catalog CHARM 
(Richichi et al. , [2005) or an estimated value taken from 
the lists of stars of reference published by Borde et al. 
(2002) or Merand et al. I|2005(l . These catalogs are indeed 
included in SearchCal. But in the general case no object 
in these catalogs complies with the specific requests of 
ASPRO in coordinates and magnitude range, and then 
the angular diameter is usually not known. A surface- 
brightness versus color-index relation should be used to 
compute angular diameter from photometry. Such rela- 
tions exist in the literature (di Benedetto 119981 van Belle 
1999J Kervella et al. I2004f) but are determined only for a 
particular luminosity class or only for few a photometric 
indices. 
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Fig. 2. Difference between measured colors of the stars 
in the Ducati catalog ( 20O and that computed for the 
same spectral type and luminosity class. The rms of this 
O-C are given for three ranges of spectral type. Only the 
dwarfs are plotted in this figure, as an example. 



Our goal was then to obtain a universal relation work- 
ing for the all luminosity classes. For that purpose we 
used, on one hand, linear diameters and an absolute mag- 
nitude determined in eclipsing binaries (which are in gen- 
eral dwarfs) and, on the other, angular diameters mea- 
sured from interferometry or lunar occultation (generally 
for giants) and apparent magnitude. 

The angular diameter 9 of a star of linear diameter D* 
(in unit of solar diameter D@) at the distance d (in pc) is 
given by: 



= Cst 



D 



(7) 



where Cst = 9.306 mas corresponds to the angular diam- 
eter of the sun seen at 1 pc. The distance of the star is 
usually a function of the apparent my and absolute My 
magnitudes: 



Then, we define the quantity ipy as: 
£>» 9 



10 (5-M v )/5 9 _306. io-™v/5' 



0) 



where ipv is computed as a function of 9 and my for stars 
with angular diameter measured from interferometry or 
lunar occultation and as function of and M v for eclips- 
ing binary components. Then, ^y-versus-color-index rela- 
tions were determined for the whole index of BVRIJHK 
system using a polynomial fit for each color index (CI). 



ipv = T, k a k CI k 



(10) 



(8) 



To determine our relations we compiled (i) the stel- 
lar diameter (from interferometric measurements, lunar 
occultation and eclipsing binaries) from Barnes et al. 
(HHZSI), Andersen ijTMTj) . Segransan et al. PTOljl . and 
Mozurkewich et al. I|2003j> and (ii) BVRIJHK photome- 
try (from Ducati 2002 and Gezari et al. 119991 catalogs) for 
a large sample of stars of spectral type O to M and for the 
whole luminosity class. We plan to regularly add new and 
more accurate measurements in our compilation, and to 
refresh our relation in SearchCal. Such relations are also 
very useful for other studies and will be described in detail 
in a forthcoming paper. 

In the top of Fig.[3|we plot one example of the relation 
for the (V — K) color. The angular diameters could then 
be computed by using Eqs. 9 and 10. The (O-C) (dif- 
ference between the measured and the computed angular 
diameters) are calculated for the stars of Segransan et al. 
l|2T7n3"jl and Mozurkewich et al. The distribution of 

the relative (O-C) are shown in the bottom of the Fig. [3] 

The first cause of uncertainty in the computation of the 
angular diameter is the variance of the calibration resid- 
uals (variance of the (O-C)). This variance includes the 
intrinsic dispersion of the relations, the error on the mea- 
sured diameters and the error on the photometry. The 
three relations with the best accuracy are ipy(B — V), 
ipy{V — i?), and ipviV — K) with an uncertainty of 8%, 
10%, and 7%, respectively. We choose these three relations 
to determine the stellar angular diameter in SearchCal. 
The polynomial fit of these three relations are given in 
Table H 

The error in the photometry is propagated in the an- 
gular diameter when a surface brightness relation is used. 
As already mentioned, we impose the restriction that the 
K magnitude of the calibrators are from observations and 
that the stars are in the Kharchenko (2001) catalog, where 
the B and V magnitudes are present. Thus, for the three 
colors used to determine the angular diameter, only the 
(V — R) color is determined from our spectral type - color 
relation, then could have substantial error. In Fig. 01 we 
show the propagation of the errors in (V — R) on the an- 
gular diameter. 

In conclusion, our angular diameter determination has 
an accuracy of <10% if measured photometric data are 
used. When the angular diameter is computed with a (V — 
R) calculated data, the accuracy is ~20%. 
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col. Ind. 


Validity domain 


ao 


ai 


a2 


&3 


db4 3,5 


Accuracy 


B-V 


[-0.4; 1.3] 


0.33822617 


0.76172888 


0.16990933 


-0.0803159 


0.36842746 — 


8% 


V-R 


[-0.25; 2.8] 


0.29974514 


0.90469909 


-0.0438167 


2.32526422 


-1.4324917 0.43618476 


10% 


V-K 


[-1.1; 7.0] 


0.32561925 


0.31467316 


0.09401181 


-0.0187446 


0.00818989 — 


7% 



Table 2. Polynomial coefficient of the ipv = T,kdkCI k relation for the three color-index retained. This relation is only 
defined for a given validity domain in color. 
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Fig. 3. top: ipv as a function of (V — K) color. The red 
line is our polynomial fit. bottom: Distribution of the an- 
gular diameters (O-C) from the ipv versus (V — K) colors 
relation. The red curve is a Gaussian function fitting the 
distribution of the (O-C). A9 is the difference between an- 
gular diameters computed from our relation and measured 
angular diameters 



For each star, a coherence test of the photometry is 
done by comparing the computed diameters with the dif- 
ferent color indexes, 6[BV], 6[VR], 9[VK]. The star is 
rejected from the list if one value of the angular diameter 
differs from more than 2a from the mean value. 

In the top of Fig. 3, the enhanced scattering of stars 
with (V — K) > 4 (spectral types later than M0) around 
the polynomial fit mainly reflect an intrinsic dispersion 



Fig. 4. Error in the determination of the angular diameter 
when the ipv versus (V — R) relation is used and when 
the (V — R) color have errors of -0.2, -0.1, 0.1, and 0.2 
magnitude. 



of the photometric data for cool evolved stars. This can 
be related to stellar variability or to the presence of a 
circumstellar envelope and a color dependent variation of 
the computed angular diameter results. Then this type of 
star cannot be considered as a good calibrator, so it is 
rejected from the final calibrator list. 



4.4. Computation of the squared visibility 

The visibility of the calibrator must be computed rigor- 
ously to avoid any differential effect with the visibility of 
the scientific target. For wide spectral band interferome- 
try, a monochromatic visibility can be computed for an ef- 
fective wavelength taking the spectral-energy distribution 
of the star into account (calibrator or scientific target) 
across the filtered spectral bandwidth. For spectrally re- 
solved interferometric measurements (wavelength-resolved 
fringes), the polychromatic visibility must be computed 
across the observed spectral bandwidth. In this version 
of SearchCal, we give only an estimation of the visibility 
for the central wavelength of the photometric band used 
for the observation. Each star in the list is considered as a 
uniform disc. The squared visibility V^ al and its associated 
error AV^j are computed as a function of the angular di- 
ameter #(mas) and its error A9 for the given instrumental 
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configuration (wavelength A(nm) and maximum baseline 
B max (m)): 

V c l^\2Mx)/x\ 2 (11) 

Ay^,=8J 2 (s)|J 1 (a;)/x|Ae/fl (12) 

with x = 15.23B max 8/ X. 

To compute the visibility, the value of the diameter 
can either be the measured 4> u d or <pid taken from the cat- 
alog CHARM, the computed <pid given in the Catalog of 
Calibrators for Long Baseline Stellar Interferometry, or 
the photometric angular diameter (/)[VK]. Using Qu in- 
stead of 9 u d in Eq.[ll] induces a bias in the computed 
visibility, SV 2 = V ca 2 (8 ud )~V ca 2 (9id) which must be esti- 
mated. With the linear representation of the limb darken- 
ing, a good approximation of the ratio 6id/9 u d as function 
of the limb darkening coefficient u is given by: 

0M/0ud=[(l - u/3)/(l - 7u/15)} 2 (13) 

with 0.0 < u < 1.0 and then 1.0 < 9i d /0 ud < 1.12. The 
bias on the visibility is maximum for a full darkened disk 
(it = 1.0) and it is always less then 5% for x < 1.0 and 
V 2 al > 0.75 (see Table EJ. 

Table 3. x max as a function of the maximum bias 6V 2 on 
the computed visibility 







0.05 
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0.025 


0.65 


0.01 


0.4 



Then, for a requested accuracy AV 2 of the computed 
visibility, the bias SV 2 can be neglected if the angular 
diameter of the potential calibrator satisfied the condition: 

8(mas) < x max (6V 2 )X(nm)/15.23B(m) (14) 

Figure [S] shows the values of 6 (mas) as function of the 
base length B and of the maximum bias SV 2 for the K 
band (A = 2.2/im). 

5. Technical aspects 

SearchCal has been designed as a distributed application 
that: 

— retrieves user requests (either through a GUI interface 
or via an command line); 

— shifts through CDS-based stellar catalogs to retrieve a 
large number of stellar parameters, according to vari- 
ous scenarii; 

— computes the missing photometry, if necessary using 
the relations mentioned in Sect. 4; 

— classifies stars as potential calibrators; 

— presents the list to the user through a GUI interface 
and handles further requests for search refinement and 
sorting. 
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Fig. 5. The maximum value of the angular diameter al- 
lowed for a calibrator as a function of the baseline for 
different values of the expected visibility bias at X2.2fim. 

SearchCal reads and writes star catalog data in VO- 
Table format (standardized XML-based format defined for 
the exchange of tabular data in the context of the Virtual 
Observatory and return format of the CDS requests). 
The resulting VO-Tables are parsed to select the possible 
calibration stars using the "libgdome" XML parser. Its 
Graphical User Interface (GUI) is a Java applet running 
on the client side in any web browser. It has been devel- 
oped using XML to Java Toolkit and is fully integrated 
in the JMMC's ASPRO Web software, allowing the user 
to display, sort, filter and save the catalog of calibration 
stars. The server side application is written in CH — h us- 
ing a flexible and scalable object-oriented methodology. 
The design allows easily the application to be updated to 
follow the improvements in scientific knowledge (in, e.g., 
the scenarii), as well as changes in the web queries and 
evolutions in data formats. 

The input panel of searchCal is presented in Fig. 7. 
As already described, the parameters of the request are: 
the observing wavelength, the range of magnitude of cal- 
ibrators in the observing photometric band, and the field 
(maximum distance of calibrators in right ascension and 
declination). The output panel is presented in Fig. 8. In 
the "result" window, a summary of the output of the re- 
quest is given as three numbers: the number of stars re- 
turned from CDS request as potential calibrators, the ini- 
tial number of potential calibrators selected after the co- 
herence test of the photometry, and finally the proposed 
calibrators without variability or multiplicity indications 
in the Hipparcos catalog. The final list of calibrators is 
displayed in the central window. The origin of each pa- 
rameter is encoded by a color code corresponding to the 
name of the relevant catalog. For the computed parame- 
ters (missing magnitude, angular diameter, visibility), the 
color code indicates the confidence level based on the qual- 
ity of the photometric data. The detailed description of the 
tables, as well as the function of the different buttons, are 
given in the SearchCal Help available as a PDF file in the 
ASPRO web site. 

Finally, from the final list of potential calibrators, one 
can refine the selection of the calibrators by changing, a 
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Fig. 6. Input panel of SearchCal 
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Fig. 7. Output panel of SearchCal 



posteriori, the parameters of the request: field around the 
scientific target, object - calibrator magnitude difference, 
spectral type and luminosity class, accuracy on the cali- 
brator visibility, indication of variability, or multiplicity. 

6. Conclusion 

We have described the principles of our calibrator's se- 
lection tool dedicated for optical interferometric observa- 
tions. Based on an online CDS request and a dedicated 
computing program, we built a powerful piece of soft- 
ware that is already open to the astronomical commu- 
nity. Our application (SearchCal) and the concurrent one 
(GetCalWeb) are the only software able to find suitable 
calibration stars in the vicinity of bright scientific tar- 
gets and that are based on a dynamical approach using a 
Web-based interface. However, differences of strategy and 
method can be noted between these two softwares. In the 
current version of SearchCal, we impose the knowledge of 
the magnitudes V and K, whereas in GetCalWeb, the mag- 
nitude K is deduced from magnitude V and the spectral 
type. The angular diameters are calculated in SearchCal 
by a surface brightness method whereas in GetCalWeb, 
they are calculated as black bodies of the selected spec- 
tral type and of magnitude V. 

In SearchCal, the limiting magnitude attainable for the 
selected calibrators is imposed by the magnitude of the 
fainter stars for which the maximum number of astromet- 



ric and spectro-photomctric parameters are available in 
the catalogs used for this selection, i.e. typically V mag- 
nitude < 10 or K magnitude < 5. In practice, these limits 
agree with the sensitivity of the interferometers currently 
in operation in the visible or the near infrared. With the 
gain in sensitivity expected with the instruments AMBER 
and PRIMA on the VLTI, it will be necessary to find 
fainter calibrators. We have already started the develop- 
ment of an extended version of SearchCal for K magnitude 
> 5. 

The use of UCD and VOtable, as well as the split- 
ting of our API in three modules ( "Access to CDS" , 
"Computation", "Display"), will allow us to easily con- 
tinue the development in the framework of the Virtual 
Observatory concept. Development of a CDS web service 
and display in the environment of a VO portal are foreseen 
in near future. 
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Appendix A: Spectral type - luminosity class - 
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Table A.l. Adopted colors for the dwarfs, in Johnson 
system. 
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Table A. 2. Adopted colors for the giants, in Johnson sys- 
tem. 
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.45 


.13 





.36 


.06 





GO 


.76 


.85 


.51 


.19 





.40 


.07 





G2 


.87 


.99 


.58 


.22 





.47 


.08 





G4 


.97 


1.07 


.65 


.28 





.50 


.08 





G5 


1.02 


1.11 


.67 


.32 





.53 


.09 





G6 


1.06 


1.12 


.67 


.32 





.53 


.10 





G8 


1.15 


1.15 


.69 


.30 





.54 


.11 





KO 


1.24 


1.24 


.76 


.34 





.58 


.12 





Kl 


1.30 


1.32 


.80 


.35 





.62 


.13 





K2 


1.35 


1.40 


.86 


.37 





.66 


.14 





K3 


1.46 


1.58 


.94 


.42 





.72 


.14 





K4 


1.53 


1.77 


1.04 


.48 





.76 


.15 





K5 


1.60 


2.10 


1.21 


.61 





.99 


.15 




K7 


1.63 


2.14 


1.22 


.63 





.99 


.16 





MO 


1.63 


2.17 


1.23 


.65 





.97 


.17 





Ml 


1.63 


2.27 


1.27 


.63 


— 


1.02 


.17 


— - 


M2 


1.64 


2.44 


1.33 


.64 




1.03 


.18 




M3 


1.64 


2.79 


1.47 


.72 




1.07 


.19 




M4 


1.64 


3.39 


1.73 


.86 




1.02 


.20 




M5 


1.62 


4.14 


2.17 


.89 




1.02 


.25 





Table A. 3. Adopted colors for the supergiants, in 
Johnson system. 



